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I. INTRODUCTION 

Solid state realization of q-bits that do not decohere 
easily is a challenging task in the field of Quantum 
Computation. Topological defects in strongly correlated 
quantum many body systems are protected from deco- 
herence and have been suggested as q-bit candidatesi"— . 
In this context, Kitaev constructed a remarkable two 
dimensional quantum spin model that exhibits abelian 
and non-abelian anyons and is exactly solvable for its 
spectrum.®'^. It was later shown that the spin correla- 
tion functions are also exactly solvableS.. This model is 
also extremely interesting form point of view of frustrated 
spin models and the physics of the resonating valence 
bond (RVB) states^rJ^. It reahzes, in an exact fashion, 
the phenomenon of quantum number fractionization and 
emergent gauge fields that were conjectured and approx- 
imately realized in 2D models for RVB states or quan- 
tum spin liquidsi^ii^. It has also been shown that the 
Jordan- Wigner transformation in this model yields a lo- 
cal fermionic theory^^"— . This makes the Kitaev model 
an important one that warrants further investigation; 
and no wonder that an extensive body of research^Sr^^i 
has already been carried out exploring its many fascinat- 
ing aspects. Kitaev showed that the model has a natural 
formulation in terms of a Majorana fermion interacting 
with Z2 gauge fields. The remarkable feature of the Ki- 
taev model is that the gauge fields turn out to be static. 
This greatly simplifies the dynamics leading to the exact 
computation of the spectrum and spin-spin correlation 
functions. 

The gauge theory of spin-i models has a recent history. 
It was initially formulated™ in the context of strongly 
correlated electronic systems such as a spin- i Mott insu- 
lator, as a way of implementing the single electron occu- 
pancy constraint, 

^(cU:.-l)|^)=0. (1) 



This equation is the Gauss law constraint for a L/(l) lat- 
tice gauge theory. It resulted in a strongly interacting 
U(l) gauge theory formalism of spin-i models. Soon it 
was realized that in spin-i systems, the U{1) gauge in- 
variance always implied an SU{2) invariance^^"— . The 
consequent SU (2) gauge theory formalism was found to 
be useful in the context of relating apparently different 
mean field solutions of the model. The extended Hilbert 
space of the spin-i system is much smaller than that of 
normal lattice gauge theory where the gauge degrees of 
freedom on the links are SU{2) group elements. Conse- 
quently, it was shown that the essential physics of the 
spin-i system is captured by a Z2 gauge theor y^^'^" , 
where the Z2 gauge group is the center of the original 
SU{2) gauge group. The Z2 gauge theory formalism has 
been effectively used to bring out the physics of quantum 
number fractionization in spin-i systemal^. 

In this paper we follow the route charted out above 
in the context of the Kitaev model and show that the 
Z2 gauge theory can indeed be thought of as the cen- 
ter of the SU (2) gauge theory of RVB theory. This sheds 
light on Kitaev's assertion^ that the model represents the 
same universality class of topological order as RVB. We 
examine the degeneracy of states in the system defined 
with periodic boundary conditions in detail. We show 
that this degeneracy, which characterises the topologi- 
cal order in the system, arises from the so called large 
gauge transformations. Namely, gauge field configura- 
tions which correspond to the same fiux configuration 
that are not related to each other by local gauge trans- 
formations. These topologically distinct gauge field con- 
figurations can be labelled by the value of the Wilson 
loops that wind around the torus in the two different 
directions. We find these gauge configurations for the 
ground state and demonstrate the four-fold degeneracy 
of the ground state by explicity computing the energies. 
Further, we generalise this proof to all eigenstates by 
constructing the operators that generate the large gauge 
transformations and showing that they do not change the 
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energy of the system in the thermodynamic hmit. 

The rest of the paper is organised as folows. Section 
reviews the Kitaev honeycomb model with its 
main features and some mathematical notions which 
are used in later sections. The SU{2) gauge symmetry 
of the model is reviewed in section where we 

show that Kitaev's choice of the representation of the 
spin operators in terms of Majorana fermions amounts 
to a SU{2) gauge fixing procedure which fixes the 
gauge upto the center Z2 gauge transformations. In 
section (lIVp . we construct a generalized Jordan- Wigner 
transformation and show that it is a Z2 gauge fixing 
transformation in the Kitaev formalism. We derive 
the fermionic Hamiltonian and simplify the conserved 
quantities in terms of gauge invariant Jordan- Wigner 
operators. In section (fV)) . we explicitly work out the 
four fold degeneracies of the Kitaev model on a torus. 
We derive the fermionic spectrum for flux free sector 
in (jV Ap and show that it leads to four degenerate 
ground states in thermodynamic limit. Following this in 
section (jVBp . we demonstrate that every eigenstate of 
the Kitaev model has four fold degeneracy. To this end 
we derive four mutually anticommuting operators that 
commute with the hamiltonian in the thermodynamic 
limit. The minimum dimension of the representation 
of the 4-dimensional Clifford algebra is 4. Thus we are 
able to prove that all eigenstates are at least four fold 
degenerate. We summarize our results in the section 

(ED. 



II. THE KITAEV MODEL 

A. The Hamiltonian 

The Kitaev model is a spin-i system on a honeycomb 
lattice. The Hamiltonian is 

H ~ Jyll^^o - J^T.^t'^h (2) 

{ij)^ {io)y (ij)^ 

where i,j run over the sites of the honeycomb lattice, 
{ij)a, o- — x,y, z denotes the nearest neighbor links ori- 
ented in the a'th direction as shown in Fig. [1] We will 
be working with periodic boundary conditions which are 
defined as follows. The honeycomb lattice is a triangular 
lattice with a basis of two sites. The sites of the triangu- 
lar lattice are given by, 

Rm,n ^ mei + ne2, (3) 

where m, n are integers and ei • 62 = — ^, ei • ei = 1 = 
62 -62- The label i of the sites of the honeycomb lattice 
therefore stands for (m, n, a) where a = a,b is the sub 
lattice label. The periodic boundary conditions are then 
defined by, 

'^ni,n,a ~ '^m+M,n+N,a- (4) 




FIG. 1: Honeycomb lattice and the Kitaev model 

B. The conserved quantities 

There is a conserved quantity associated with every 
plaquette of the lattice. If the plaquette is denoted as p 
and its vertices labelled as shown in Fig. [U then, follow- 
ing Kitaev's notation, the conserved quantity is, 

B.^alalalalalal. (5) 

We have — \ implying that the BpS can take values 
±1. It is clear that any product of the BpS will also com- 
mute with the Hamiltonian. In fact there is a conserved 
quantity associated with every closed self avoiding loop, 
C, on the lattice defined the following way. At every site, 
the path will pass through two of the three bonds that 
emanate from it. We call these two bonds as the tan- 
gential bonds and the third one the normal bond. We 
associate two tangential vectors at each site, tii and t2i 
which are either x, y or z according to the direction of 
the incoming bond and the outgoing bond respectively. 
We then define a normal vector hi as 

n., EEtHXt2j. (6) 

If the sites of C are 11,12, in, then the conserved 

quantity associated with it is, 

N 

Bc=X{{n^^■a^^). (7) 

n=l 

It can be checked that, 

[Sc,H]-0, Bl^l. (8) 

We will call C topologically trivial if it can be written 
as a product of BpH. On the torus, we have two loops 
which wind the torus around in the two directions which 
cannot be expressed as a product of BpS. One cannot be 
obtained from the other by multiplication by BpS. We 
will call these two (Wilson) loops Wi and W2- 

All the BpS are not independent due to the identity, 

X{Bp = l. (9) 

p 

Thus there are Np — 1 independent BpS, where Np = MN 
is the number of plaquettes. Together with Wx and W2, 
we have a total of Np + 1 conserved quantities on the 
torus. These two loop operators account for the four fold 
degeneracy on a torus. 
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III. FROM 5(7(2) TO Za 
A. U{1) and SU{2) gauge symmetry 

Interacting quantum spin systems often lead to spon- 
taneously broken symmetric states such as a ferro, an- 
tiferro or spiral magnetic states. Low energy physics of 
these ordered states are captured by the well known spin 
wave approximations, reasulting in Goldstone mode type 
bosonic low energy effective theories'^'^. A quantum spin 
liquid, on the other hand, has no classical long range or- 
der. Experiments in LaCu04 and other low spin Mott in- 
sulators, according to Anderson, indicated possible pres- 
ence of neutral fcrmionic excitations in a quantum spin 
liquid^"—. A theory to describe such a quantum spin 
liquid or resonating valence bond (RVB) state needed a 
paradigm shift from spin wave theory. It was also clear 
that a quantum spin liquid, in view of different possible 
phase coherence among disordered spin configurations, 
could offer a variety of quantum spin liquid states to be 
realized in nature. RVB gauge theory attempted to cap- 
ture these new possibilties, through an approach involv- 
ing enlarged Hilbert space and emergent gauge fields in 
strongly correlated electron systems. Through the work 
of Wen"^^ and others it has become clear that there is a 
plethora of spin liquid phases, characterized by quantum 
order and projective symmetry groups. 

The spin-i Hilbert space can be realized as the sub- 
space of the Hilbert space of two fermions defined by the 
constraint, 



=n 



(10) 



where cj. and Ca are the fermion creation and annihi- 
lation operators. As mentioned earlier, Eq. ([TOl) can be 
looked upon as the Gauss law constraint for a U{1) gauge 
theory. The LHS of the equation being the generator of 
the following U{1) gauge transformations on the fermion 
operators, 
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The spin operators, 



(11) 



(12) 



are then the gauge invariant observables of the theory. 
ct" are the Pauli spin matrices. 

The single occupancy constraint in the spin-i theory 
implies that a spin-f hole is the same as a spin- J, particle 
in the physical space. This can be mathematically ex- 
pressed as an SU{2) gauge invariance. It is convenient 
to express this symmetry in terms of a matrix of the 
fermion operators 



(13) 



In terms of this matrix, the spin operators are given by. 



5« = itr 
4 



(14) 



The generators of the SU{2) gauge transformation are 
given by. 



S" = -itr ^a''^. 
4 



(15) 



The ^ matrix transforms under the SU{2) spin and 
SU{2) gauge transformations as. 



(16) 



Where 17s and Ug are SU{2) matrices representing the 
spin and gauge transformations respectively. It is clear 
from equations ((Ti)) . and (|16p that the spin oper- 
ators are gauge invariant and the generators of gauge 
transformations are spin singlet. The constraint in 
Eq. pU)) is exactly equivalent to the SU{2) Gauss law. 



s^ii^) =0. 



(17) 



Before we close this section, we wish to mention that 
the above gauge theory formalism offers a possible way 
to understand quantum spin liquid states as and when 
they exist. This formalism does not gaurantee a simple 
gauge theory structure at all energy scales in the physics 
of the problem. It only suggests that in some systems 
(for some Hamiltonians) at low energy scales there could 
be emergent gauge fields and interesting consequences of 
quantum number fractionization, quantum order etc. At 
high energy scales gauge fields interact and it is no more 
simple or useful to talk in terms of emergent gauge fields. 
As we will see soon, the Hamiltonian invented by Kitaev 
on a honeycomb lattice is very special. It offers static Z2 
gauge fields and makes the Z2 gauge theory meaningful 
at all energy scales. 



B. Majorana fermions and the Z2 theory 

We can make connection to Kitaev's representation of 
the spins by writing. 



Ct = n ' = 



(18) 



where c, Cx, Cy and Cz are Majorana fermions. The single 
occupancy constraint, Eq. (ITUl) reduces to exactly Ki- 
taev's form, 



CCxCyCz — 1. 



(19) 



Kitaev's representation of the spins then get written as. 



-CCa = S^~S^ 



(20) 



Note that these three operators are not equal to the 
gauge invariant spin operators in the extended Hilbert 
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space but are exactly equivalent to them in the physical 
Hilbert space. Substituting the expressions in Eq. (^0)) 
for the spin operators in the Hamiltonian is then equiva- 
lent to adding gauge fixing terms. Since the SU (2) gauge 
generators are invariant under the Z2 center of the gauge 
group, these terms only fix the gauge upto the central 
Z2 group represented by e^^^^ . The Hamiltonian will 
therefore continue to have a Z2 gauge symmetry. 

The simple example of a spin-i in a magnetic field 
illustrates these issues. If we take the Hamiltonian to be 

H = B S^. (21) 

then the theory has SU{2) gauge symmetry, the degen- 
erate ground states in the extended space are 

\GS) = acjlO) +/3|0) +744IO), (22) 

for arbitrary a,/3 and 7, |ap -I- + I7P = 1. The last 
two states in the RHS of Eq. ((22|) transform as a dou- 
blet under the SU(2) gauge symmetry. Gauge averaging 
therefore projects out the ground state in the physical 
subspace, namely cj^|0). 

If the Hamiltonian is taken to be, 

H ^iB cc^^B {S^ - S^^ . (23) 

The theory has only Z2 gauge invariance, the degenerate 
ground state in the extended Hilbert space are, 

\GS) =ac\\{)) +-ic\c\\Q). (24) 

The second state in the above equation transforms non- 
trivially under the Z2 gauge transformation. Thus again, 
under gauge averaging, the ground state in the physical 
sector is projected out. Thus Kitaev's representation of 
the spin operators can be interpreted as adding gauge 
fixing terms to the SU{2) gauge invariant hamiltonian 
which leave an unbroken (unfixed) Z2 gauge symmetry. 

IV. THE Z2 GAUGE THEORY OF THE KITAEV 
MODEL 



The link variables are then given in terms of the occu- 
pancy number of the bond fermions, 

'^(u)a = 2 x{,j)„X(»j% - 1- (27) 
It is easy to see that, 

= 0. (28) 

Thus the link variables can be interpreted as static Z2 
gauge fields. 

It is remarkable that at one shot Kitaev hamiltonian 
has been solved exactly for the entire many body spec- 
trum ! Infact, two related phenomena occur: i) the Z2 
gauge theory is exact at all energy scales and ii) the en- 
larged Hilbert space gets decomposed into sectors that 
are identical gauge copies having the same energy eigen 
values (Fig. [2]). The Hilbert space enlargement does not 
produce any unphysical state, but only gauge copies. In 
the standard U(l) RVB gauge theory for a Heisenberg 
antiferromagnet, for example, it is easy to see how un- 
physical states are brought in by Hilbert space enlarge- 
ment. For example, an unphysical state containing M 
doubly occupied and M empty sites gives spectrum of a 
Heisenberg antiferromagnet containing 2M missing sites. 
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" Z Gauge Copes 





Number of lattice sites: 2N 
Dimention of Physical Hilbert space: 2 

2N 

Dimension of the Enlarged Hilbert space: 4 

Kitaev Hamiltonian has identical eign value spectrum in each Box 

FIG. 2: How the enlarged Hilbert space for Kitaev Model 
become just Gauge Copies. 



A. The Hamiltonian 



Following Kitaev^, we write the hamiltonian in terms 
of the Majorana fermions. 



H 



2^ 2^ *C,W(y)„Cj, 

1=1 {ij)a 



(25) 



where the link variables are defined as = iCaiCaj- It 

is natural to express them in terms of the bond fermions^ 
defined as, 



(26) 



B. The Jordan- Wigner transformation and Z2 
gauge fixing 

A remarkable feature of the Kitaev model is that the 
Jordan- Wigner transformation yields a local fermionic 
Hamiltoniaoi^"— . In this section we show that the 
Jordan- Wigner transformation in the Kitaev model can 
be interpreted as a Z2 gauge fixing procedure result- 
ing in gauge invariant (gauge fixed) Majorana and bond 
fermions. The choice of the Jordan Wigner path amounts 
to a choice of the Z2 gauge. 
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1. The Jordan-Wigner fermionisation 



The gauge fixed Hamiltonian 



We define the Jordan-Wigner transformation as foi- 
lows^. Take any Hamilton patli on tire lattice defined 
by a sequence of sites i„, n = 1, . . . , Ng, where Ns is 
the number of sites in the lattice. The path will classify 
each bond as normal or tangential as defined in section 
IIIBI The normal bonds will form a dimer covering of 
the lattice. For a given site we attach three vectors, 
two tangential vectors denoted by tii, t2i and one normal 
vectors fii, such that they follow Eq. ([6]). Then we can de- 
fine the Jordon-Wigner transformations in the following 
compact way, 



2. Gauge invariant Jordan- Wigner fermions 

To define the Jordan-Wigner transformations we as- 
sociate two Majorana fermions (77^,^ and ^i^) at a given 
site 'i„'. These Jordan-Wigner (Majorana) fermions are 
defined in terms of the gauge invariant spin-operators in 
the following way , 



n(fi. 



(29) 
(30) 

(31) 



m — 1 



It can be easily checked that the above definitions refer 
to the usual anti-commutations relations for Majorana 
fermions, 



(32) 



Since the Jordan-Wigner fermions are constructed en- 
tirely from the spin operators, they are manifestly gauge 
invariant. However, it is also interesting to see this by 
rewriting equations (l30l) and (|29|) in terms of the original 
Majorana fermions and gauge fields. 



iCi 



1) ilii ■ Cii- (33) 
■ -Ui^h) iih ■ Cii. (34) 



The transformation can be inverted to write the spins 
in terms of the fermions. 



(Ji. 
CTi. 



(35) 

(36) 
(37) 



The Hamiltonian can be written in terms of the gauge 
invariant fermions as, 

(«i>x («i>B (y% 

(38) 

where the gauge fixed Z2 fields, Uij are, 

Uij — normal bonds, 

= 1 tangential bonds except (ij) = (iiijvs), 



Ns 
n=l 



Ns 
n=l 



iS is a well known conserved quantity in the one- 
dimensional applications where it corresponds to the to- 
tal number of fermions modulo 2 and determines the 
boundary conditions on the fermions. Thus the Jordan- 
Wigner transformation is equivalent to a gauge fixing 
procedure where all the gauge fields on the tangential 
bonds (except one) are set equal to 1. The choice of the 
Hamilton path amounts to a gauge choice since it defines 
which of the bonds are tangential. It also defines the sign 
in the definition of Uij for the normal bonds. In Eq. p9p . 
the sign corresponds to a Hamilton path which winds 
regularly in the ei direction as shown in Fig ([3]). Here- 
inafter all explicit computations will be with respect to 
this path. A general algorithm to go to this gauge, which 
we will refer to as the Jordan-Wigner gauge, is given in 
appendix El 




FIG. 3: Jordan-Wigner path on a torus for 4x4 lattice. The 
numerics at the lattice sites describe how the Jordan-Wigner 
path traverses the lattice. 



4- The fermionic conserved quantities 



These completes the definitions of Jordan-Wigner 
fermionisations used in this article. 



All the gauge fixed Z2 fields are conserved quantities. 
It is convenient to express them in terms of the gauge 
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fixed bond fermions on the normal bonds, 



V. FOUR FOLD DEGENERACY ON A TORUS 



^ ^ i^ + i ^ ^^ - ^^3 



(40) 



The conserved quantities are then the occupation num- 
bers of the bond fermions, 



2 XijXij - 1- 



(41) 



There are Np normal bonds and hence the bond fermion 
occupation numbers form a set of Np conserved quan- 
tities. Thus along with S, we have Np + 1 conserved 
quantities consistent with the analysis in terms of spin 
variables in section III Bl To see the meaning of 5, it 
is convenient to define complex fermions in the matter 
sector from the two ^ fermions on every normal bond. 



irij 



,t = ^» - 
Vij — 



2 ' 2 
It can then be shown that, 

^0 = E ^l-^^j-' 



(42) 



z bonds 

E 

normal bonds 



(43) 



5. From Kitaev gauge to Jordan-Wigner gauge 

We have explained that Jordan-Wigner gauge is a spe- 
cial realisation of Kitaev gauge where all the gauge fields 
residing on the tangential bonds are fixed to unity. One 
may wonder whether there exists a gauge transformation 
on the lattice which renders Eq. ([25]) to that of Eq. (p8| . 
Indeed there exists such a gauge transformations. Refer- 
ring to the Jordan-Wigner path given in Fig. we do 
the following gauge transformations at a site 'n' 



n 



i+ic„, 1 < n < iV. 



(44) 



In the above equation the various indices corresponds 
to the Jordan-Wigner path in the Fig. Ui^j+iS are 

the Z2 gauge fields that would normally exists on the 
link joining site and 'i + 1' if we apply the Majorana 
fermionisation as adopted by Kitaev. After implementing 
the above gauge tranformations Z2 gauge fields appear 
only on the normal bonds. The above gauge transforma- 
tions yield a new conserve quantity S' which appear on 
the bond where the Jordan-Wigner end path meets. The 
expression for 5' is, 



5' 



-1 



(45) 



= 1,N 



The significance of S' is very similar to S introduced in 
Eq. (123) • 



Since the gauge fields are static, the problem is re- 
duced to one of non-interacting fermions on a lattice. It is 
knownSiSI that the lowest fermionic ground state energy 
is obtained for the flux free configuration, i.e. Bp = 1, \/p. 
On the torus, there are four gauge in equivalent configu- 
rations for every configurations of Bp as argued in IIIBI 
These correspond to the four values of the gauge invariant 
conserved quantities Wi and W2 ■ First We examine these 
four different gauge field configurations for fiux free sector 
and compute the corresponding fermionic ground state. 
Following this we show that this leads to four fold de- 
generacy of ground states in thermodynamic limit. Next 
we demonstrate explicitly how to obtain the four fold de- 
generacy for every configurations of fluxes. To this end 
we derive the required operators which enables us to ob- 
tain any one of the inequivalent gauge field configurations 
from the other for arbitrary flux conflgurations. 



A. Degenerate ground states on a Torus 

To start with we briefly recapitulate the no- 
tions of Jordan-Wigner transformation and refer to Fig. 
([3]). The normal bonds are the ones that form the ba- 
sis of the triangular lattice except for the (0,n, a) line. 
On this line the normal bonds are the ones between 
(0, n,a) and (0,n + l,b)). We choose the first site of 
the path to be ii = (0,0,6). The four flux free conflg- 
urations are then explained as follows. To this end we 
write the exact Hamiltonian and Bp for this particular 
realisation of Jordan-Wigner transformation. We divide 
the Hamiltonian in three parts. Hint , Hbound and H^nd- 
Hint includes all the internal bonds and Hbound includes 
all the boundary bonds except one where the Jordan- 
Wigner end points meets. Hend includes the interaction 
for the bond where the Jordan-Wigner end points meet 
each other. Similarly all the Bp's are categorized in the 
above three different way. Below we write the various 
parts of the Hamiltonian and Bp's. 



Hint — E] ^'^a:^m,n'7m+l,n-(-l + "^JyVn 



(46) 



^ — ' m,n+l 
m,n 



Where u m,n is defined on each boundary y-bond. The 

m,n+l 

Hamiltonian for the end bond is given by, 



(48) 
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Now with the definition of ■0 fermion and x fermion we 
get the equations (|47p to be rewritten as, 

H^nt = ^ Ja;(^^,„+^m,n)(V'm+l^„+l -V'm+l,n+l) 
m.n 

+ Jz (2^1„,„^„,„ - 1). (49) 
The Hamiltonian for the slanting bonds, 

Hbound = y^Jy"^ ('/'m,n+'/'m,n)(V'm,„+l -■0m,n+l) 

^ — ^ m.n+1 

m , 71 

+ J.{2^l^^^m,n-l). (50) 

Lastly the Hamiltonian term for the end bond where the 
end points of the Jordan- Wigner path meet each other is 
given by, 

Hend = -SJz{2%l;\,j j^^M.N - I)- (51) 

Following equations ([M)) . (15(11) and P^ . we rewrite the 
complete expressions for various conserved quantities ap- 
peared in the final form of the Hamiltonian. The conserve 
quantity Um,n defined on each internal z-link is given by, 

Um,n = (2xL,«Xm," " !)• (52) 

Similarly the conserved quantity defined on each bound- 
ary y-bonds (which are labeled by the z-bonds it is con- 
nected with (i,e u m,n )) is given by, 

m,n+l 

U m,n = (2x^m,n X ^1)- (53) 

And the S, for the Jordan- Wigner gauge, is given by, 

5^ _(_1)MJV+N,,+W,^ (54) 

If and denote the parity operators for the gauge 
fermions and the matter fermions respectively, then we 
can write S = -{-lY^'^V^V'^ . From the Fig. H we 
see that a single hexagon always contains two normal 
bonds where each normal bond is associated with a con- 
served(static) Z2 gauge field . Bp for any plaquette is 
the product of these two conserved Z2 gauge fields. Thus, 
Bp = UijUki , where 'ij' and '/cZ' are the normal bonds for 
the plaquette 'p'. For the end plaquette where Jordan- 
Wigner path terminates Bp is given by Bp — —Suijiiki- 
Now we are in a position to show the ground state de- 
generacy in thermodynamic limit. To this end we ex- 
plicitly write the four in-equivalent gauge field configura- 
tions corresponding the flux free configurations and write 
down the corresponding fermionic Hamiltonian. Finally 
we find the spectra for each of these four fermionic Hamil- 
tonian. 



1. Choice 1 

Here the flux free conflguration is obtained by mak- 
ing all the u's to be 1 and S = —1. The loop conserve 
quantities are having the following eigenvalue ,Wi = 1 
and W2 ~ 1 . These particular choice makes the resulting 
fermionic Hamiltonian transitionally invariant and usual 
Periodic boundary condition in both the direction can be 
used to diagonalize the Hamiltonian. We will explicitly 
write the Hamiltonian in terms of complex fermion i/;. 
However one can equivalently work in terms of 7y, Ma- 
jorana fermion representation. To keep this in mind we 
will continue to mention appropriate gauge transforma- 
tions for T] fermions as well as t\} fermions. The complete 
translational invariant Hamiltonian is given by, 

B = Xl^^(^m,"+^™'»)(^m-(-l,n+l -VAn+l.n+l) 
m.n 

+ -^2/((V'm,n + V'm,n))(^I,^„+i " V'm,n+l) 

+ J,(2V'1„^„Vm.n-l). (55) 

This is a manifestly p-wave superconducting Hamil- 
tonian and can be easily diagonalized by going to the 
momentum space. The constraint on the number of ?/; 
fermion becomes, 

n(2V'^^„V'm.n - 1) = 1, (56) 

which implies that we are to fill up only the even number 
of fermions. We define the Fourier transform of the 
i>m,n as given below, 

where k\ — ^ttjj and k2 — ^n-^. This is obtained by 

noticing the fact that we can write k — jjGi + j^G2 

where G1/2 are the reciprocal lattice vectors which are 
given by, 

^ 47r / V3 1 \ ^ 47r 
G,^—\^—e^ + -e,j ; = ^e,. (58) 

Substituting this we get the resulting Hamiltonian in mo- 
mentum space as, 

k 

+ J2^k-MNJ,. (59) 
fc 

Where ek = 2(Ja;C0sfc^ -I- Jy cos ky + J^) and 5k = 
2( sin fcj, 4- Jj^ sin fcj,). k^ — k.rij,, ky — k.rij, and n.x^y = 

ie^. ± ^Gy are unit vectors along x and y type bonds. 

In Eq. (1591) . the sum over 'fc' runs over first half of 
the Brillouin zone and does not include the 'fc'-points 



8 



(tt, 0), (0, tt), (0, 0). The first line of tlie Hamiltonian is 
diagonalize by the following transformations, 

( ak \ _ ( cos6'k -i sin^k \ ( V'fe \ 

V /3fe y ^ -i sin^k cos(?k ) \ VAfe ) • ^ ' 



Where cos26Ia; f-k/Ek, with Ek = + 5^.. Then 
re- witting the Hamiltonian we get, 

H = J2Ek{alak-l3ll3k) 

-^(eo.o + eo.TT + Ctt.o) + ^efc' - iV^ J3)(61) 

fc' 

Here sum over k' runs over full Brillouin zone. The last 
term in the parenthesis is always zero for Torus. 

S. Choice 2 

Here the flux free configuration is obtained by making 
all m's -1 and S = —1. The corresponding values of loop 
conserve quantity is given by,VFi = — 1 and W2 = — 1. 
To implement the Fourier transformation we do the fol- 
lowing steps. We make the following gauge transforma- 
tion. r]\,[ ^ — —f]Mn for all 'n'. In terms of ip fermion 

the necessary gauge transformation is 'ipAi^n — ^''/'l/ for 
all 'n'. Then the Hamiltonian requires r]M+i,n = —Vi.n 
and rim,N+i = Vm.iior alternatively 4>M+i,n = -i'l^n 
and iprn,N+i = ipm,!)- This is equivalent to anti-periodic 
boundary Condition in ei direction and periodic bound- 
ary condition in 62 direction. The necessary Fourier 
transform is defined with, 

with fci = ||(m -I- i); k2 — ^n^. Substituting this in 
the Hamiltonian and diagonalizing straightforwardly we 
get, 

k k' 

Here also 'A:' runs over first half of the Brillouin zone 
and the 'fc" runs over the full Brillouin zone. Note the 
absence of (0, tt), (tt, 0) (0, 0) mode. They do not appear 
here for this anti-periodic boundary condition. This will 
be true for the choices 3 and 4 also. Various parameters 
appearing in Eq. (|63p are given below. 

El = ^Jiel+Sl), 

ek = 2{Jx cosfca; + Jy cosky - Jz), 

6k = 2{Jx sin kx + Jy sin ky). (64) 



3. Choice 3 

For this case the fiux free configuration is obtained by 
making iiM.n = ^1 where 'n' runs from 1 to A'^ — 1. All 
other {t's are 1 and 5 = 1. The loop conserve quan- 
tity Wi takes value 1 and and W2 -1. Similar to the 
previous case we need to do following gauge transfor- 
mations in order to apply Fourier transform. We make 
'?m N = N f^'' tcrms of fermion the nec- 

essary gauge transformation is tpm,N = n fo'' all 'to. 
The resulting Hamiltonian requires rim,N+i = ^Vm,i and 

riM+l,n = ril,n{ i>m,N+l = -^m,l and V'M+l,n = i'l^n)- 

This indicates anti-periodic boundary condition in 62 di- 
rection and periodic boundary condition in ei direction. 
The resulting Fourier transform is defined with, 

m , 2t: , 1 , , ^ 

A:i=27r— ; fcs = — (n+-). (65) 

The resulting Hamiltonian in 'k' space is similar to 
Eq. ([63|) with = 2{Jx cos k^ + Jy cos ky + Jz). 



4- Choice 4 

For this choices we need, UM,n — 1 where 'n' runs 
from 1 to iV — 1 . All other it's are -1 and 5 = 1. Wi = 
— 1 and W2 = 1. In this case one requires combined 
gauge transformation mentioned for the choices 2 and 
3. This makes the Hamiltonian anti periodic in both 
the directions. The required Fourier transform is defined 
with, 

27r, 1, , 27r, 1, 
^1 = ^(^+2); *2 = ^(ri + -). (66) 

Proceeding as before we get exactly Eq. (|63|) with iden- 
tical expression for efc. 



5. Ground State Energy in thermodynamic limit 

To get the ground state energy for the gauge choices 2,3 
and 4, one fills up the negative energy states consistent 
with the boundary condition (i.e, to satisfy the constraint 
S which restricts the total number of particles (odd or 
even number) to be taken). In the limit of M,N — 00 
all the above four choices, the ground state energy is 
obtained as, 

Eg^^I E{k,,ky)dkidk2, (67) 
167r^ Jbz 

with E{kx, ky) is defined before. The appearance of a '-' 
sign in the expression of ek for choices 2 and 4 can be 
accounted for by shifting the ki integral to tt — /ci . Thus 
it is clear that in thermodynamic limit ground state has 
four fold degeneracy. 
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B. Four fold Topological degeneracy for any 
eigenstate 

In the preceding section, we have shown that four de- 
generate ground states are characterized by four different 
topofogicaUy distinct gauge field configurations. From 
the section IIVI we infer that every flux conflguration is 
characterized by such four topofogicaUy distinct gauge 
field configurations. This leads to a four-fold degeneracy 
for every eigenstate, including the ground state, in the 
thermodynamic limit. 

We now demonstrate this explicitly. Our proof is simi- 
lar to that of Wen and Niu"^^ which shows the topological 
degeneracy of fractional quantum Hall states on a torus. 
We construct two operators which we call Vi and V2, 
that act on states with a given gauge field configuration 
and produce states with a different gauge field config- 
uration without changing the values of the flux opera- 
tors. Bp. They however change the values of the Wilson 
loop operators, Wi^2- These two operators are therefore 
the generators of the so called "large gauge transforma- 
tions". There are four topofogicaUy different sectors of 
gauge field configurations corresponding to = ±1. We 
further show that only effect of the large gauge transfor- 
mations on the matter sector is to change the bound- 
ary conditions of the Majorana fcrmions from periodic 
to anti-periodic or vice versa. Thus the energy eigenval- 
ues only change by ~ where L is the length of the 
torus. The eigenstates in the four sectors (related by the 
action of Vi) are therefore degenerate in the thermody- 
namic limit. 

The four operators Vi, Wi characterise the topological 
degeneracy. We show that Vi and Wi satisfy the following 
algebra, 

{Vi,V2} = 0, {V,W,} = 0, [V,Wj],^j = 0. (68) 
We can then construct four operators, Tq, a — I, ... ,4, 

Ti^ViWi, T2 = V2W2, T3 = Vi, Ti^V2. (69) 
These satisfy the Clifford algebra, 

{Ta,Tb} = 25ab (70) 

Thus we show that the four-fold topological degeneracy 
on the torus is characterised by the 4-dimensional Clifford 
algebra. 

We will first write down and discuss the expressions 
for Vi and V2 for the 32 site system illustrated in Fig. 
|3l The construction easily generalises for any even-even 
lattice. 

^1 = at<jl<jla^r (71) 

7 23 

V2 = olol.ol.ol^'^^al^ol (72) 

i=2 1=18 



The two Wilson loops for this lattice are, 

M/ — ^y ^y ^y ^y ('Ti\ 

8 

W2 = n< (74) 

1=1 

It can be verified that the above constructions satisfy the 
algebra in equation and hence the topological oper- 
ators defined in equation ((69)) satisfy the Clifford algebra 
(I7ni) . It can also be verified that Vi commute with all the 
Bp's. 

Now consider simultaneous eigenstates of the Wilson 
loop operators, 

Wi\wi,W2) = Wt\wi,W2) , Wi ^ ±1 (75) 

The algebra in equation (j68|) implies, 

Vi\wi,W2) ^ \ ~ Wi,W2), V2\wi,W2) ^ \wi,-W2) (76) 

Thus we have shown that Vi are the generators of large 
gauge transformations. 

We next consider their action on the hamiltonian. 

ViHV{'^ = H\ V2HV2^ = (77) 

where is the same as H, except that the bonds on 
one non-trivial loop in the ei direction have changed 
sign, namely the loop (1, 2, 3, 4, 5, 6, 7, 8, 1). In H^, a line 
of parallel bonds in the ei -I- 62 direction have changed 
sign. Namely, the bonds (7 - 8), (31 - 32), (23 - 24) 
and (16 — 17). We will now write down the operators 
for a general even-even lattice and then show that the 
transformed hamiltonians are degenerate in the thermo- 
dynamic limit. We will show that in the fermionised the- 
ory, these changes of sign can be absorbed into the sin- 
gle particle eigenfunctions of the Majorana fermions and 
change the single particle energy eigenvalues by ~ 1/L. 
Thus the energies do not change in the thermodynamic 
limit, making every many body eigenstate four-fold de- 
generate. 

The general expressions for T^i(2) are given by 

M 

^1 = n ^m!o' 
m=l 

N/2 m=l,M-l 

n=l n=l,N/2 

It can be verified that these constructions satisfy the alge- 
bra in equation (j68[) and also commute with all the Bp's. 
In , the bonds (m, 0, b) — (m -f 1, 1, a) and (m, 1, a) — 
(to, 0, h) change sign and in , the bonds (Af — 1, n, 6) — 
(0, n + l,a) change sign. 

In the fermionised theory, the single particle eigenfunc- 
tions satisfy the equation, 

Y,zA,,^-^e-d,^ (79) 
3 
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where Aij is an antisymmetric matrix coupling the near- 
est neighbours of the honeycomb lattice. The eigenvalues 
come in pairs and we denote (</>")* — e^" — — e". n 
will then go from 1, . . . , NM. The Hamiltonian is diago- 
nal in terms of the complex fermions defined by, 

= + (80) 

i 

= ^e"(2ata„-l) (81) 

n 

We now make the transformation, 

C™,a = e-'^/^C™.. (82) 
Equation ([79]) then gets written as, 

J2 = E i'^'j + = ^"'^r' (83) 

j 3 

where is the antisymmetric matrix corresponding to 
and 5 Aij oc 1/N when N is very large. Thus the 
single particle energy eigenvalues of A and A^ are iden- 
tical in the thermodynamic limit when TV — cx). The 
spectrum of H and are also therefore identical. The 
mapping of the eigenvalues of H and can be similarly 
shown using the transformation, 

= e™-/*V:^™,a (84) 

We can also write equation ((78)) in terms of the gauge 
invariant fermions. Vi or V2 can be written as Vi — 
V/^Vf {i = 1,2). As the notation suggests, V/^'s and 
Vf 's are composed of operators which belong to the mat- 
ter sector and the gauge sector respectively. They are 
given by, 

VF^l[Cn,o, m=l,M-l, 

N/2 j = l.N/2 

n=l k=2,N 
N/2 m=l,M-l 
^2^ — W ^0,2n^0,2n W '7m,2n-l^m,2n-2-(85) 
n=l n=2,N/2 

From the above expressions we can easily find the com- 
mutation relations of parity operator for gauge fermions 
and the matter fermions. We find that, 

{7'^,yi}-0, [7'^,^2] = 0, [V^,Vi] = Q, [V^,V2]^Q. 

(86) 

The fact that the parity of the matter fermions are 
conserved is consistent with the fourfold degeneracy dis- 
cussed here and in (|V A|) . 



VI. DISCUSSION 

In this paper we have discussed many important as- 
pects of Kitaev model. We have shown how the SU2 
gauge contained in the Kitaev model and explained ex- 
plicitly how this SU2 gauge symmetry is reduced to Z2 
gauge symmetry. We have solved the Kitaev model us- 
ing Jordan- Wigner method in a general way. Though 
Jordan- Wigner transformations has been used earlier to 
solve for Kitaev model, our formalism reveals many new 
features, for example, our definition of J-W transforma- 
tion is applied to a torus. We focused on the gauge field 
contents of the Kitaev model and the issue of topological 
degeneracy. We showed that ground state is four fold de- 
generate on the torus in both phases. While it indicates 
non-trivial topological order, the ground state degener- 
acy does not distinguish between the gapless and gapped 
phases. Finally we have shown the equivalence between 
the fermionised Hamiltonian obtained in Kitaev gauge 
and J-W gauge. Lastly we have constructed four mu- 
tually anti-commuting operators on a torus to illustrate 
explicitely the four fold degeneracy for every eigenstate. 
Our analysis reveals that Jordan- Wigner analysis can be 
used even in quantum spin liquid problems, as a gen- 
eral method, to bring out non trivial gauge field content, 
thereby providing an alternative method in resonating 
valence bond theories. 
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Appendix A: The gauge fixing algorithm 

In this appendix we detail the gauge fixing algorithm 
to go from any gauge field configuration to one in the 
Jordan- Wigner gauge. First we will write Hamiltonian 
as obtained from the Eq. ([55]) . We notice that for to ^ 0, 

And for to = 0, 

Then starting with the Ferromagnetic Hamiltonian we 
get, for TO = 0, n = 0, 

- Oao = -^vlovUS- (A3) 
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For TO 7^ 

- 0-°"^ a'"' = -in'' n'' fif" f'' ) (A4) 
now for m — 0,n ^ 

-<:n<,u^<.A..n- (A5) 

The above three equations gives complete description 
of all the z-z interaction. Now for the y-bond we get for 
TO = 

For TO 7^ we get 

At last we write Hamiltonian for x-interaction. For to ~ 


And for to 7^ we get, 

Now we make the following gauge transformation 
which we call Jordan- Wigner gauge. 

VL.n ^ i-^rvL.n- (AlO) 

Then we choose for each normal bond Uij = — i^f^j. 
This gives the Jordan- Wigner Hamiltonian given by 



equations [46l[47| and [48l Then we define complex fermion 
Xm,n on each normal internal z-bond in the following way, 

). (All) 

Similarly on each normal slanted y-bond which is joined 
with a z-link (to, n) and (to, n + I) we define, 

Cn-l-l = (X ™>» +X^m,n ) ; = — (X -X^m,n )■ 

m,n+l m,n+l I m,n+l m.n+1 

(A12) 

With this we have alway Ui,j = {2x1 jXi-j ~ At last 
on each z-link we define ijj in the following way, 

C« = ('/'m^n+V'm,™) ; Vrn,n = j (V'm.n " V"™,™) ■ (Al3) 

Now the quantity S for lattice of dimension (M, N) , 
is given by, S = --(-_i)A^JV+7V^+7V^^ jj^^.^ ^^^^ 

are the number of ip and x fermions respectively. Now 
noticing the number of gauge transformations needed 
for various flux free configuration we find that we need 
to fill even number of ■0 fermions for each different 
gauge choices. All the results derived here are based 
on a representative lattice of dimensions M and N 
where M and N are both even. It is straight forward to 
carry the analogous calculation for lattice where M and 
N can be anything, odd or even. However the results 
obtained here should not change in thermodynamic limit. 



^ M. A. Nielsen and I. L. Chang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, England 2000). . 

^ C. H. Bennett and D. P. DiVinccnzo, Natme 404 
247(2000). 

^ A. Yu. Kitaev, A. H. Shen and M. N. Vyalyi, Classical and 
Quantum Computation (America! Mathematical Society, 
2002). 

* S. Das Sarma, M. Friedman and C. Nayak, Phys. Today, 

59 32(2006). 
^ J. Preskill, Phys. Today, 52 24(1999). 
^ A.Yu. Kitaev, Annals of Physics, 321, (2006), 2-111. 

A.Yu.Kitaev, Annals of Physics, 303, (2003), 2-30. 
® G. Baskaran, Saptarshi Mandal, R Shankar, Phys. Rev. 

Lett, 98, 247201(2007). 
^ X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 

11413(1989). 

° N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773(1991). 
^ D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 
2376(1988). 

^ N. Read and B. Chakraborty, Phys. Rev. B 40, 7133(1989). 
^ R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 
1881(2001). 

^ V. Kalmayer and R. B. Laughlin, Phys. Rev. Lett. 59, 



2095(1987). 

Senthil. T, Fisher. M.P.A, Phys. Rev. B 62, 7850(2000). 
Xiao-Yong Feng, Guang-Ming Zhang, Tao Xiang, Phys. 
Rev. Lett. 98, 087204(2007). 

Han-Dong Chen and Jiangping Hu, Phys. Rev. B 76, 
193101(2007). 

Han-Dong Chen and Zohar Nussinov, J. Phys. A: Math. 
Theor. 41, 075001, 2008. 

Saptarshi Mandal, R Shankar, G Baskaran, Poster 14, In- 
ternational Conference on Mott Insulator, IISc, Bangalore, 
July,2006, http://www.physics.iisc.ernet.in/seminars-past.phpj 
K. P. Schmidt, S Dusuel and J. Vidal, Phys. Rev. Lett. 
100, 057208(2008). 

K. Sengupta, D. Sen, and S. Mondal, Phys. Rev. Lett. 100, 
077204(2008). 

G. Baskaran, D. Sen and R. Shankar, Phys. Rev. B 78, 
115116(2008). 

S. Mandal, N Surendran, Phys. Rev. B 79, 024426(2009). 
S. Dusuel, K. P. Schmidt, J. Vidal, Phys. Rev. Lett. 100, 
177204(2008). 

G. Baskaran, Z. Zou and P. W. Anderson, Solid St. 
Commn. 63973(1987) 

Ian Affleck, Z. Zou, T. Hsu and P.W. Anderson. Phys. Rev. 
B 38 745(1988). 



12 



E. Dagotto, E. Pradkin, and A. Moreo. Phys. Rev. B 38 

2926(1988). 

G. Baskaran. Indian J. Phys. 80 (6), 583-592(2006). 

G. Baskaran, E. Tosatti and L. Yu, Int. J. Mod. Phys. B, 

2, 555, 5(1988). 

J. B. Martson, Phys. Rev. Lett. 61, 1914(1988). 

E. Pradkin, Field Theories of Condensed Matter Systems, 

HarperCollins Canada / Perseus Books. 

E. Licb, F. Wu, Phys. Rev. Lett. 20, 1444(1968). 

V. Kalmeyer, R. B. Laughlin, Phys. Rev. Lett. 59, 

2095(1987). 



P. W. Anderson, Science 235, 1196(1987). 
P. W. Anderson, Science 288, 480 (2000). 
X. G. Wen, Phys. Rev. B 65, 165113(2008). 

E. H. Lieb, Phys. Rev. Lett. 73, 2158 (1994). 

X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377(1990). 
Abhinav Saket, S. R. Hassan and R. Shankar, Phys. Rev. 
B 82, 174409 (2010). 

F. J. Burnell and Chetan Nayak, Phys. Rev. B 84, 125125 
(2011). 



